COEXISTENCE OF THE CRITICAL SLOWING DOWN AND GLASSY 
FREEZING IN RELAXOR FERROELECTRICS 
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We have developed a dynamical model for the dielectric response in relaxor ferroelectrics which 
explicitly takes into account the coexistence of the critical slowing down and glassy freezing. The 
application of the model to the experiment in PbMgi/3Nb2/303 (PMN) allowed for the reconstruc- 
tion of the nonequilibrium spin glass state order parameter and its comparison with the results 
of recent NMR experiment ( Blinc et al., Phys. Rev. Lett 83, No. 2 (1999)). It is shown that 
the degree of the local freezing is rather small even at temperatures where the field-cooled permit- 
tivity exceeds the frequency dependent permittivity by an order of magnitude. This observation 
indicates the significant role of the critical slowing down (accompanying the glass freezing) in the 
system dynamics. Also the theory predicts an important interrelationship between the frequency 
dependent permittivity and the zero-field-cooled permittivity, which proved to be consistent with 
the experiment in PMN (A. Levstik et. al., Phys. Rev. B 57 11204 (1998)). 



I. INTRODUCTION 

Relaxor ferroelectrics which are disordered per- 
ovskite ferroelectrics like PbMgi/3Nb2/303 (PMN) or 
PbSci/2Tai/203 (PST) represent a new class of materi- 
als which have been a subject of numerous investigations 
(see, for example, Rcfs. |l|- |l^] ) . Relaxor ferroelectrics 
manifest themselves in the extraordinary low frequency 
dispersion of their dielectric permittivity compared with 
regular ferroelectrics. The position and the height of 
the permittivity maximum plotted as a function of tem- 
perature depends on the frequency of the probe field 
and shifts to lower temperatures when the frequency de- 
creases. 

The latter behavior is accompanied by the observed 
splitting between the field-cooled (FC) and zero-field- 
cooled (ZFC) permittivity and the existence of long lived 
remanent polarization. However, as has been emphasized 
recently, the bulk of the relaxation spectra in relaxors re- 
mains active even far below the temperature Tf where the 
FC and ZFC permittivity split, and the nonlinear suscep- 
tibility does not diverge at Tf. All these findings indi- 
cate nonequilibrium phenomena and quasi-nonergodicity 
( "freezing" ) , rather than a true thermodynamic dipolc 
spin glass transition. 

On the other hand, the existence of the very high di- 
electric constant indicates that these systems are close to 
ferroelectric instability, and, therefore, one could expect 
the manifestation of the critical slowing down of dynam- 
ics and a competition between the critical slowing down 
and the dipole spin glass freezing. Such a dichotomy 
makes it a non-trivial task to extract from the experi- 
ment quantitative characteristics of the freezing. 

One of the important characteristics of spin glass freez- 
ing is a value of the Edwards- Anderson order parameter 
q, which in relaxors is a time dependent quantity. We 



will show that this quantity can be extracted from the ex- 
perimental results on the frequency dependent permittiv- 
ity, with the use of the model discussed below. Also the 
theory predicts an important interrelationship between 
the frequency dependent permittivity and the zero-field- 
cooled permittivity, which proved to be consistent with 
the experiment. 



II. MODEL 

It has been recently proposed |l2[ that relaxor behav- 
ior is a common characteristic of the collective dynamics 
of localized giant dipole moments distributed in highly 
polarizable crystals, and a necessary condition of relaxor 
behavior is the simultaneous existence of the broad dis- 
tribution of the local field and the broad distribution of 
dipole relaxation frequencies. Also experiments |10| , pd] ,p[ 
favor the hypothesis that the physical origin of localized 
dipole moments is the off-center shift of the atom even 
at high temperatures. In this model the large values of 
the dipole moments are associated with the polarization 
cloud (cluster) formed by the simultaneous displacements 
of the other atoms adjacent to a given off-center ion. 

A convenient approach to describe dynamical behavior 
of relaxor ferroelectrics is to start from the Bloch type 
equations widely explored in the theory of regular ferro- 
electrics possessing Debye relaxation. We write the Bloch 
equation in the form 



dPci 
dt 
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Eq.(|l]) describes the relaxation of polarization of each 
cluster to its quasiequilibrium value P^i{E) which de- 
pends on the value of the local field E induced by other 
clusters at any moment of time. In general the local field 
_E is a time dependent random field. It includes also the 



contribution from the applied field Eex and the contri- 
bution from the static random field caused by material 
imperfections. Note that Eq.(|l|) is quite general and, al- 
though the explicit form of P^l^ {E) and the precise defini- 
tion of the polar clusters are model dependent, it effects 
only the coefficients of the theory discussed below. 

In order to apply Eq.p) to relaxor ferroelectrics one 
should perform an average over the distribution of re- 
laxation times T and the random local field E. We as- 
sume that the distribution function f{E, P) of the local 
field, which depends parametrically on the average po- 
larization of the system P{t), has the form f{E,P) = 
f{E - 7P - jeoEex/'iTT). This form of f{E) is con- 
sistent with the mean field approximation /(i?, P) = 
S{E — jP — 7eoi?e2:/47r), where 6 is the Dirac S- function 
and 7 is the local field phenomenological parameter. The 
effect of the local field fluctuations corresponds to the re- 
placement of the (5-function by the function / with finite 
width. The value jeoEex/^iTr is the local field induced by 
the external field at the location of each off-center ion in 
a dielectric media with the dielectric constant eo >> 1, 
i.e., we assume that the localized dipole moments are dis- 
tributed in crystals with large lattice permittivity caused 
by the existence of the soft modes. This assumption is 
consistent with the observation in relaxors of the soft 
optical modes (remaining finite at all temperatures) , fre- 
quencies of which are of the same order of magnitude (~ 
SOcm"-'^) as those in highly polarizable dielectrics such as 
KTaOa, SrTiOa. 

It is known that the existence of the broad distribution 
of relaxation times leads to non-exponential behavior in 
the polarization relaxation. In order to reproduce this 
effect within the proposed formalism we rewrite Eq. M) in 
the integral form and then take the average with respect 
to r, E, and the initial cluster polarization Pci{0). Thus, 
we obtain 



Pit) = P(0)g(i) - k{T) 1^ dt'^[P{t t') 



dt' 
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In Eq. (0) we assumed a linear response of the polarization 
to the applied electric field, and T > T^ where T^ is the 
temperature of a possible ferroelectric phase transition. 



k{T)= f dEP- 



e,^j,^ df{E,P) 



dP 
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and for second order phase transitions A; — > 1 as T —> Tc, 
The function q{t) is equal to 



q{t) = / dT~g{T)e-"^ 



drgir) 



(4) 



where ^(t) is a distribution function of relaxation times 
and the right hand side expression of Eq. (0) is valid for 
smooth functions ^(t). The variable q{t) describes the 



fraction of clusters effectively frozen at time t and there- 
fore has the meaning of the Edwards- Anderson spin glass 
order parameter on a finite time scale. Note that, as we 
assumed above, q{t) — > at i ^ cxd. 



III. FREQUENCY DEPENDENT PERMITTIVITY 

The steady state frequency dependent permittivity 
can be easily obtained from Eq.(H) assuming Eex{t) — 
-Eex e*"* and using the definition 



e(w,r) =47r 

Thus we obtain 

e(u;,r) = 

where 

Q(w,T) 



aPH 



eo 



eo- 



\-k{T)Q(uj,T) 



dre" 



dt 



(5) 



(6) 



(7) 



For the smooth function ^(t) the real part oiQ{uj, T) can 
be simplified as |14 



Q'{cu,T)^ 



dt- 
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The relation(g) between the functions q{t) and Q'{uj) will 
be employed below to obtained information on the degree 
of local freezing in relaxors from the frequency dependent 
dielectric constant using the fact that usually Q" < Q' . 
At ijj = Eq.(^ defines the static or field-cooled per- 
mittivity 



^FC 



eo 



1 - kiT) 



(9) 



On the other hand, at LUTm >> 1 the function Q{lu, T) — > 
and therefore e{uj) —> eq- Thus, in our model eo can be 
regarded as a high frequency permittivity e^o ■ Combining 
equations (0) and (ph we obtain 



Qiu) 



e(c^) ^ ep 
^FC — eo 



e_Fc 
eo 



(10) 



Note that in a number of recent publications on re- 
laxors |l^,|l^,Q the analysis of the relaxation spectrum 
Q{uj) was based on a different relation between Q{ll>) and 
the permittivity, namely 



Qiu) 



e(c^) - ep 
eFC ~ ep 



(11) 



where we substituted ep for eoo- Equation (nil) was in- 
troduced earlier in spin glasses, assuming that the relax- 
ation time T in Eq. (H) is a characteristic of independent 



cluster-relaxators. However this equation has a different 
meaning when appHed to relaxor ferroelectrics. Indeed, 
Eq. ( pdl ) can be obtained from Eq. (^ if we first perform 
the average over the distribution of the local random 
fields at a constant value of t and then perform the av- 
erage over T. With such a two-step averaging we arrive 
at Eq. (pi]) where r is replaced by t* = re/eo- Note that 
r* is a relaxation time of the collective polar mode with 
the wave vector (7 — » 0, undergoing the critical slowing 
down of dynamics, rather than the relaxation time of the 
individual dipoles or clusters. Thus the applicability of 
Eqs.(pl|), and (||) to a system with long range or meso- 
scopic polar order implies that the crystal can be divided 
into macroregions within which the relaxation time of all 
dipole moments has the same magnitude, and the aver- 
age over T means the average of the dielectric response 
for different macroregions. On the other hand, Eq.jlQ) 
is consistent with the formation of the short range clus- 
ters in the absence of correlations between the relaxation 
times of different clusters. 



IV. ZERO-FIELD-COOLED SUSCEPTIBILITY 

AND REMANENT POLARIZATION 

The values of ezpcit) can be obtained by solving 
Eq.(p|) with P(0) = and E^x — const. We consider here 
a particular case when one can approximately neglect the 
memory effects in Eq.(2) by replacing P{t — t') « P{t) ( 
which implies a fast decay of dq/dt). Thus, we obtain 



ezFcit) = — 



eo 



Ml - Q{t)) 



(12) 



One can see that ezFc{t) is identical to e(w) given by 
Eq.(g), if one substitutes w hy 1/t in Eq.(||). 

In order obtain the remanent polarization Pr(t) we 
solve Eq.(H) in the same manner as above, but with 
E^x = and P(0) 7^ 0. We obtain 



m 



P{t) 

P(0) 1 - k{l - q{t)) 



^FC — ^ZFC 
^FC -^ £0 



(13) 



The denominator in Eq. (|l3| ) reproduces the effect of the 
critical slowing down of dynamics which is imposed on 
the effect of glass freezing characterized by the slow de- 
caying function q{t). For example, when /c — > 1 ( i.e., 
the system is in the vicinity of the second order phase 
transition temperature ) it follows from Eq.(|l3|) that 
Pr{t) — const independently of the value of q{t). 



COMPARISON WITH THE EXPERIMENT IN 

PMN 



A. Nonequilibrium spin glass order parameter 

The nonequilibrium spin glass order parameter q{t^ T) 
is a very important quantity determining the dielectric 
response of relaxors. The dielectric permittivity can be 
formulated in terms of the parameter g(i, T) in a way 
consistent with the description of magnetic susceptibility 
in magnetic alloys where spin glass and ferromagnetic 
order coexist. It is convenient to rewrite Eq.(0) in the 
identical form 



e(u;,r) 



fceo(l - q) 
l-k{l-q) 



eo, 



(14) 



which separates the contribution to the permittivity from 
the localized dipole moments and the lattice permittivity 
eo- The first term in Eq.([l4|) is very similar to the well 
known solution for the susceptibility given by the infinite 
range Sherrington-Kirkpatrick model (mean field theory) 



X = 



C{l-q) 
T-9{l-q) 



(15) 



Equation (^_5|) has been widely used for the description 
of the experiments on ZFC susceptibility as well as fre- 
quency dependent susceptibility in magnetic alloys flj] 
and relaxors ||l^ by treating the parameters C and 9 as 
purely phcnomcnological fitting parameters. In fact, our 
approach justifies and generalizes Eq.(^ in the case of 
nonequilibrium spin glass state. The deviation of the pa- 
rameter k(T) from 9 /T indicates the deviation from the 
mean field picture. 

Equations (||) and (g|) allow reconstructing the val- 
ues of the functions q{t, T) from experiments on the fre- 
quency dependent permittivity. We applied them to the 
recent data on PMN by Levstic et al |13| who obtained 
the values of the frequency dependent permittivity, as 
well as field-cooled and zero-field-cooled permittivity. 
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FIG. 1. Reconstructed values of the nonequilibrium spin 
glass order parameter in PMN as a function of temperature 
corresponding to the following frequencies of the applied field: 
lOOkHz (1), IkHz (2), and 20Hz (3). Squares reproduces the 
data obtained in Ref.[19]. 

We used eo ~ 1000 corresponding to the values of e 
in the low temperature limit. The results of the recon- 
struction are presented in Fig.l where the temperature 
dependences of q{T) as shown for three distinct times 
corresponding to the frequencies of the appHed field used 
in the experiment equal to 20Hz, IkHz, and lOOkHz. One 
can see that the values of the q{T) are rather small even 
at T=200K where the values of e((jj) are almost 10 times 
less than the values of epc- The explanation of this im- 
portant effect lies in the extremely high dielectric con- 
stant of PMN leading to fc « 0.98 at T =200K. Thus, 
the system is extremely close to the ferroelectric insta- 
bility which results in the significant effect of the critical 
slowing down, accompanying the glass freezing, upon the 
system dynamics. 

In Fig.l the reconstructed values of q are compared 
with the the data obtained from the recent experiments 
by Blinc at al. |19|] on the observation of inhomogeneous 
broadening of NMR lines of Nb ions. The appearance 
of inhomogeneous broadening signals the slowing down 
of the dynamics on the scale of inhomogeneous linewidth 
(equal to 10-50 kHz). A reasonable correspondence be- 
tween curve 1 and the data |l^ (with the account of 
comparable time scales in both experiments) allows us 
to suggest the following scenario for the NMR broaden- 
ing mechanism in PMN. As known, e.g., from the anal- 
ysis of nuclear spin lattice relaxation of Li and Nb ions 
in KTaOarLi, Nb |Q, the appearance of polarization on 
the mesoscopic length scale in highly polarizable crystals 
results in a significant modulation of the electric field gra- 
dient which leads to the modulation of nuclear resonance 
frequency. Since in disordered PMN crystals all atoms 
lack a center of symmetry, one can assume the following 
form of the nuclear frequency shift of Nb atoms 



Ai^ = aPr 



(16) 



where lS.v is the frequency shift of the atoms belong to the 
given polar cluster. The width 5 of the inhomogeneous 
line will be proportional to the fraction of polar clusters 
( given by the nonequilibrium spin glass parameter q(T)) 
whose dymamics is effectively frozen on the time scale of 
1/(5. The above broadening mechanism gives also a nat- 
ural explanation of the Gaussian line shape observed in 
Ref. [|l9[ since the distribution function of cluster polar- 
ization is close to Gaussian in highly polarizable crystals 
pO| , and, according to Eq. ( [iq ) , the shape of inhomoge- 
neous NMR line coincides with the shape of the distribu- 
tion function of cluster polarization. 



B. Interrelationship between ZFC and frequency 
dependent permittivity 

Another practical aspect of the proposed model, which 
helps to clarify the experimental situation in PMN, is the 
interrelationship between the frequency dependent per- 
mittivity and zero-field-cooled permittivity obtained in 
Sec. 4. In order to employ this interrelation note first that 
in PMN, at the temperatures nearby the temperature of 
the permittivity maximum, the characteristic relaxation 
time satisfies the Vogel-Fulcher (FV) law P,H . It allows 
to replace r in Eq.(Q) by 



— T„pT-Tn 



Toe 



(17) 



and introduce the temperature independent distribution 
function g{U) of potential barriers |ljjlj]. We obtain 

{t > To) 



q{t) = 1 - 



(T-To)ln(^) 



dUg{U) 



(18) 



Thus, <7 is a function of (T — Tq) ln(t/To) that leads to the 
validity of the scaling relation q{t, T) = q{ti, Ti) with 



Ti = (T - To 



H-) 
H-) 



-To 



(19) 



Using Eqs.(p^, (g), and (^) one can reconstruct the val- 
ues of q{t,T) and Q'{lj,T) and, therefore, the values of 
e(w) by employing the experimental data on (.(lui) at a 
given frequency wi. 

The reconstructed values of e(ci;) at a; = 2QHz are pre- 
sented in Fig. 2 (solid line) where the data at wi = IkHz 
have been taken as a reference. The parameters of the 
fit are Tq « 223K and tq ~ 4 • 10~^^sec, which are al- 
most the same as obtained in Ref. [Of To ~ 224K and 
To « 4.3 • 10~^^sec. One can see that at high tempera- 
tures {T > 24:5K) the reconstruction reproduces the ex- 
perimental data with good accuracy. The deviation be- 
tween the reconstructed and experimental data for lower 
temperatures is due to violation of the VF law in PMN 
already discussed in Refs. g,|l^. The self-consistency of 
the model has been tested also by reconstructing the val- 
ues of the zero-field cooled permittivity using Eq. ( [l2| ) and 
the same reference data as above. In order to find the 
time t entering Eq. (|l2|) we use the fact that in the experi- 
ment [Q the dielectric constant was measured by slowly 
heating the crystal from 80K at the rate 0.5 K/min. It 
gives e.g., the estimate t =330 min at T = 245K. The 
result of the reconstruction is shown in Fig. 2 by the solid 
line. As one can see from Fig. 2, the position and the 
hight of the maximum of ezFC £^re reproduced with very 
good accuracy. The deviation between the reconstructed 
values and experimental data is due to violation of the 
VF law taking place at low temperatures. Indeed, the 



temperature of the maximum of ezpc corresponds to 
T w 245K, which, as discussed above, is the boundary 
value for the vahdity of the VF law. If the VF law would 
be valid for all temperatures then both ezpc and e(w) 
would approach the value eo at T = Tq. One can see 
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FIG. 2. Experimental (squares) and reconstructed (solid 
lines) values of the dielectric permittivity in PMN. Experi- 
mental values are reproduced from Ref. [13] 

that this tendency is reproduced by the interpolation of 
the solid curves in Fig. 2 to lower temperatures, thus clar- 
ifying the origin of the observed low temperature devia- 
tions between the reconstructed and experimental values 
for ezFC- 



VI. CONCLUSION 
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The high dielectric constant of relaxor ferroelectrics 
signals that these materials are close to ferroelectric in- 
stability which manifests itself in the critical slowing 
down of dynamics. In disordered relaxor materials the 
critical slowing down is accompanied by glass like freez- 
ing or cluster dynamics, which could be characterized by 
the non-equilibrium spin glass order parameter. We have 
developed a dynamical model of the dielectric response 
in relaxors which explicitely takes into account the co- 
existence of critical and cluster dynamics and allows for 
separating them from the experimental observation. 
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